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ABSTRACT 


Certain space shuttle missions may require the retrieval of 
passive spinning and processing satellites. One proposed means of 
retrieval utilizes a free-flying teleoperator launched from the shuttle. 
The feasibility of nulling the combined spin and precession (passivation) 
of a typical rigid satellite is first established using a Lagrangian “ 
formulation. It is shown that a proposed asymmetric teleoperator cannot 
be used for passivation in its present form because it would quickly 
tumble over after being spun -up to synchronize with the angular rates 
of the satellite. In addition, a dynamic analysis is made of the 
combined teleoperator-satellite system where an initial misalignment 
of their respective angular momentum vectors is assumed. 



Chapter I 
INTRODUCTION 


With the development of the space shuttle system, space tug, and 
remote manipulator units (RMU’s), it will soon be possible to maneuver 
satellites and other masses in orbit. Typically, such maneuvers are 
needed for the assembly of space stations, deployment of payloads, and 
retrieval of satellites for repair, refurbishment or removal. A 
satellite that is to be retrieved might well have some initial 
rotational angular momentum. Hence, for this satellite to be retrieved 
successfully, it may be imperative to be able to null its angular rates. 
This process is referred to as "passivation," A number of shuttle-based 
missions have already been identified which may require the retrieval of 
satellite payloads. These payloads include: 

i) Stabilized, normally -operating satellites which require periodic 
servicing or refurbishment. Stabilization of these satellites could be 
accomplished passively by spin-stabilization, or actively by reaction 
jets, control moment gyros (CMG*s), or inertia wheels. 

ii) Freely -tumbling and/or spinning satellites. The Micrcaneteorite 
Exposure Module (MEM) , for example, is currently conceived of as a non- 
stabilized passive satellite. Close inspection, removal and retrieval 
by the shuttle may be a requirement. 

iii) Satellites that have developed an attitude control system (ACS) 
malfunction. The nature of the malfunction may cause the satellite to 
be in a partially or completely unstabilized condition. 

Typical satellites to be launched in the late 1970 's and early 
1980 's for which retrieval might be required are listed in the study 
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conducted by Bell Aerospace Company. Criteria used in selecting 
satellites for study and their characteristics are tabulated in that 
study. The Research Application Module (RAM) is the roost deroanding^^^ 
in terms of size and inertia for retrieval missions. This payload is 
of maximuni possible size for the payload bay, and thus, in 5 >oses high 
accuracy requirements on shuttle positioning and stabilization of the 
teleoperator. The current configuration of the RAM indicates that it 
is to be CMC stabilized. However, gyro failure could cause motion 
about one or more axes during retrieval. 

Dynamic state of a satellite being retrieved can range from a 
completely stabilized one to the worst case of combined three -axis 
tumbling. A satellite is said to be spinning when its angular momentum 
vector is along its longitudinal axis, and is said to be tumbling when 
this vector is along an arbitrary axis, iFigure 1 shows a satellite 
coning about its angular momentum vector. The coning angle is 
influenced by the inertia properties of the satellite. Angular velocity 
can be expressed at any instant of time as the resultant of a spin 
component along the longitudinal axis and another component along the 
angular momentum vector S, as shown in Figure 2. The angular momentum 
vector also could be expressed as a resultant vector of two corresponding 
components. If one of these components is zero, the satellite is in 
either a state of symmetric spin or flat spin. Utilizing an RMU to null 
these types of spin is straightforward. The RMU approaches the satellite 
along the axis of angular motion (i.e. along the angular momentum vector), 
then docks and applies torques to remove the angular momentum. The 
problem, however, is far more complicated for the general case when the 
satellite has both spin and tumble con^onents of the angular velocity. 






4 


In this case, it is not only necessary to synchronize the teleoperator 

with the spinning and nutation rates of the satellite, but to adjust 

the grappling arm to the nutation cone angle. 

A limited amount of study has already been done to analyze "free- 

flying” and "shuttle-attached" teleoperators for retrieval purposes. 

This work has determined broad performance requirements needed by the 

teleoperators in order to retrieve typical satellites from orbit. Based 

on this, a baseline "free-flying" teleoperator concept was established 

fl 21 

and a preliminary plan developed for its implementation.'* * Specific 
objectives included determination of teleoperator performance require- 
ments for typical satellite configurations and for typical dynamic 

conditions expected during the encounter. A study conducted by Faile 
C3) 

et al . ' considered a more complicated case of dynamic passivation of 

a spinning and processing object, by applying various torques to the 
satellites. In most of these studies, it is assumed that the momentum 
vector of the teleoperator coincides with that of the satellite. This 
is, however, an oversinplification of the situation encountered during 
an actual maneuver. Firstly, it is difficult to synchronize the 
nutation and spin rates of the teleoperator arms with those of the 
satellite. Secondly, it is difficult to align the teleoperator exactly 
along the satellite momentum vector for final rendezvous. Finally, it 
is difficult to determine the exact alignment of the satellite momentum 
vector itself. Thus, it would be desirable to study the resultant 
motion of the combined teleoperator-satellite system for a certain 
degree of misalignment between the two momentum vectors. This has been 
achieved here by allowing the angular rates and the nutation angle of 
the teleoperator arms to differ slightly from the predetermined, fixed 
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values for the rates and the angle of the satellite. The misalignment 
between the two momentum vectors is then calculated and the combined 
body analyzed for its resulting motion. 

The feasibility of nulling the combined spin and precession of a 
typical rigid satellite is first established using a Lagrangian 
formulation. Next it is shown that a proposed asymmetric teleoperator 
cannot be used in its present configuration for passivation because it 
would quickly tumble over after being spun -up to synchronize with the 
angular rates of the satellite. Finally, a dynamic analysis is made of 
the combined teleoperator-satellite system where an initial misalignment 
of their respective angular momentum vectors is assumed. 
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Chapter II 
SYSTEM DYNAMICS 

The system and function requirements for each mission are discussed 
in detail in Ref. 2. System requirements include those of performance, 
information, interface, and support. The overall function requirements 
e.g., approach -rendezvous phase, capture phase, and recovery phase may 
be accomplished by either a free-flying or attached-boom teleoperator. 
Although the numerous tasks required for a free-flying teleoperator 
(FFTO) during the three phases may differ slightly in detail from those 
for an attached-boom teleoperator, they apply equally for both the cases 

2,1 Approach -Rendezvous Phase 

It is assumed that the teleoperator is carried in the shuttle pay- 
load bay during launch. A deployment mechanism is needed to provide a 
fly-off/dock station, which is positioned well away from the bay to 
assure safe operations. This mechanism could be a telescopic boom or 
a truss -type structure. When the FFTO reaches the vicinity of the 
satellite, one or two circumnavigations may be necessary in order to 
inspect the satellite and determine its angular momentum vector. This 
information is essential for aligning the FFTO in the appropriate 
approach direction for docking. 

As indicated previously, the RAM satellite was found to be the most 

demanding with respect to the function requirements. Expected yaw, 

pitch, (D and roll, u) rates for various failure modes of the RAM 
y 2 

satellite as indicated in Ref. 2 are listed below: 
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i) An abnormal shut-down of the satellite’s attitude control 
system (ACS) followed by an extended period of satellite drift would 
develop: 

(0 = 0 ) < 0.025 rad/sec. 

X y ' 

CO <0.1 rad/sec. 
z 

(ii) A failure of the ACS about one axis: 

(a) failure in roll would develop: 

(0 , (0 ^0.1 rad/sec. 

X y 

0 )^ < 10 rad/sec. 

(b) failure in pitch or yaw would develop: 

( 0 , 0 ) < I rad/sec. 

X y 

0 ) ^ 1 rad/sec. 

z 

(iii) A failure of the satellite ACS about more than one axis 
would develop: 

(0 , (0 ^ 1 rad/sec. 

X y 

0 ) <10 rad/sec. 

z 

The nutation cone angle 0 depends upon the ratio of the combined 

pitch-yaw rate (co ) to the roll rate (o) ) . This coning angle is zero 
xy z 

for a pure spin (co = 0), When the coning angle is small, the 
xy 

preferred approach for capture with convenience may be along the 
mcmientum vector to the end face of the satellite (end approach) . 
However, when this angle is large, the preferred approach may be to the 
waist of the satellite at the center of mass (waist approach) . The 
rendezvous phase ends with the FFTO about 20 to 50 feet away from the 
satellite along the final approach direction, as established by the 
inspection maneuver. 
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2.2 Satellite Motion Characteristics for the Capture Phase 

Consider a cylindrical, symmetric satellite as shown in Figure 3, 
which is both spinning and tumbling in space. When these angular rates 
are passive, i.e., there are no external torques acting on the 
satellite initially, the total angular momentum vector H of the 
satellite is fixed in direction in inertial space. The body continues 
to spin about its longitudinal axis and to process or cone about the 
angular momentum vector at a constant coning angle 9. In the absence 
of any external torques, the coning angle 0 and the spin and precession 
rates remain constant for a symmetrical satellite. 

Before one can successfully take the satellite to the shuttle, it 
is imperative to null its angular rates. This can only be done by 
eliminating the angular momentum of the satellite, i.e., nulling the 
individual components of the angular momentum vector. This is possible 
by applying torques to the body through teleoperator gripping arms in 
directions opposite to those of the w components. There are obviously 
many ways of applying torques to the body, depending upon the design 
and geometry of teleoperator arms (see Appendix for one such design). 

For a teleoperator with a freely spinning spindle at the end of the 
"hand", a passivating torque may be applied along the pitch or yaw axis 
of the satellite, which is easier than to try to reduce its spin rate. 
This observation is especially true if the teleoperator hand is designed 
to lock onto a docking ball located at the center of the satellite end 
face. Here it is obvious that the magnitude of the despinning torque 
applied by the teleoperator along the roll axis of the satellite is 
limited, as the docking ball is situated on or very near the spin-axis 
of the satellite. If the teleoperator is designed to approach the 



9 



Figure 3. Satellite Axes System. 
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spinning satellite from one end by tracking the center of the satellite 
as in the end approach , the magnitude of applied pitch or yaw torque 
would be much more than that of the rolling or despinning torque it can 
apply. The teleoperator is in a position to despin the satellite very 
effectively, if it is designed with a grappling jaw mounted on a spindle 
synchronized with the spin rate of the satellite, and the jaw grips the 
satellite at the waist. Thus, the ways to nullify components of momentum 
may be broadly classified into three approaches: 

1) eliminate either the spin component or the precession component 
and then the other, 

2) eliminate both components simultaneously by applying torques 
on spin and precession axes, or 

3) eliminate both components simultaneously by applying torques 
on spin and precession axes and a torque to simultaneously reduce the 
coning angle. This last procedure is found to passivate the satellite 
at the fastest rate. 

Suppose it is desired to eliminate the precession component of H 
first. It would then be necessary to apply a torque opposite in 
direction to that of the precession component. Application of this 
torque will obviously cause some drift of the total momentum vector, 
depending on the magnitude of the passivating torque applied. Thus, 
a torque of small magnitude will have to be applied over a long period 
of time to limit drift of the momentum vector within reasonable bounds. 

At the end of this period the precession component of fi would be 
completely removed, and the total angular momentum is then in simple 
spin. In actual practice, however, the angle 8 will not be reduced 
exactly to zero during this process. There will be some residual 
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up its arms and hand to synchronize exactly with the precession and 
spin rate of the satellite. The hand orientation is simultaneously 
adjusted to match the coning angle. The teleoperator is now in a 
position to capture the satellite smoothly, without introducing 
impact forces and torques. However, in practical rendezvous situations 
it is difficult to achieve this exact alignment of the two angular 
momentum vectors. Hence, it is also difficult to achieve perfect 
synchronization of the rates of the two vehicles. Thus, in practice, 
the teleoperator approaches the satellite along a direction slightly 
different from the direction of its momentum vector. It then tries 
to capture the satellite with its arms and hand which are spinning at 
slightly different rates than those of the satellite. It is essential 
to study the motion of the combined body after capture to ensure that 
the teleoperator does not tumble over or go into an unstable mode. 

The analysis of a system of rigid bodies, rotating at different 

( 4 ) 

angular rates and connected rigidly, is presented in detail by Martz 
and Greenwood. The solution of the second order, non-linear equations 
of motion describing the motion in terms of Eulerian angles, expressed 
in elliptic functions, is given by Whittaker. However, we derive 
the equations of motion of the composite body by Lagrangian formulation 
and solve them on a digital con5>uter by using second order Runge-Kutta 
method of integration. . 

After capturing the satellite in its arms, the teleoperator returns 
to the shuttle. This phase is known as the recovery phase. No compli- 
cated maneuvers between the two masses are anticipated. 

FACE BLANK NOT wrr.uinni 
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SYSTEM EQUATIONS OF MOTION 
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The complete equations of motion, describing the motion of the 
rigid satellite under the influence of despinning and detumbling 
torques and the motion of the combined teleoperator-satellite after 
the capture phase, are derived in this chapter. These equations 
describe the motion of the two systems in terms of Eulerian angles. 
Coordinate systems used, along with the various dynamical parameters, 
such as the components of the angular momentum vector H, the 
angular velocity vector S, etc., are shown in Figures 3 and 4. 

Figure 3 shows the coordinate system used for the satellite 
moving under the influence of passivating torques. The XYZ coordinate 
system is the moving coordinate frame in the orbit, with X-axis 
directed outwards along the radius vector) Y-axis being perpendicular 
to the orbit plane and the Z-axis lying in the orbit plane to complete 
the right handed system. The xyz system is the body fixed coordinate 
system rotating with the satellite. The z-axis is the longitudinal 
axis of the satellite along the axis of symmetry in the present case. 

The x-axis and the y-axis are any two axes, located perpendicular to the 
z-axis and passing through an origin at the center of mass. Their exact 
location is immaterial because of symmetry. The location of the body 
axes referenced to the moving frame is described with the help of three 
Euler angles. The corresponding coordinate frames for the combined 
teleoperator-satellite system is shown in Figure 4. In this case, the 
origin of the XYZ-raoving frame is located at the center of mass of the 
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. Figure 4. Coordinate Axes Frame for the Combined 
Teleoperator-Satellite System. 
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control unit of the teleoperator. The xyz-system is again the body- 
fixed system, with the x-axis (and hence the y-axis) being arbitrarily 
chosen initially. The location of the body axes is once again 
described relative to the moving XYZ-frame with the help of three 
Euler angles. 


3.1 Analysis of the Satellite Motion 

Non-constant angular velocities under the influence of the 
despinning and detumbling torques. The equations of motion of the 
satellite, moving under the influence of the applied torques, can be 
derived from the total kinetic energy T^ of the satellite via a 
Lagrangian formulation. The total kinetic energy T^ of a symmetric 
tumbling body is 


12 2 2 
T=^[I(a)^ + I (O^ + I ] 

s2 xx yy zz 

S X 

From the geometry of the Figure 3, it is seen that 


(A) = d> + Ip cos 6 
z 

2 2 2 ;2 . 2 _ ^ *2 
0 ) = u) + oj = ^ sin 0 + 6 

xy X y 


( 1 ) 

( 2 ) 

(3) 


and 


I = I 


(4) 


• '"'s = i J ’•x * I ^ 'i' 

* * . s s s 

Now, the Euler-Lagrange equations of motion for a rigid body are 

(Ref. 5) 

A SiT Of 

(6) 


d ,9L . 9L _ . 
dt ^ ‘ " 


^ 3q r 

r 


where Lagrangian L = T - U, and = generalized force in, the direction 
th 

of the r^ generalized coordinate. 
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In the present case, since we are interested in only relative 
motion of the satellite with respect to the moving frame, we neglect 
the potential energy U of the satellite. In view of this, 
differentiating the expression for the kinetic energy of the satellite 
with respect to the three generalized coordinates, i.e., the three 
Eulerian angles, will yield the equations of motion of the satellite, 
as follows: 

ip direction: 

Cl sin e + I cos'^6) + 2 (I - I ) sin 0 cos 0 6 

^s ^s 

+ I (4> cos 9-0 sin 00) - (7) 

^s ^ 

0 direction: 

® ^ ® 6 + 1 i sin 6 = (8) 

A Z X z ^ 

s s s s 

(j) direction: 

1 C<l> + cos 6 - |i0 sin 6) = t_ (9) 

s ^ 

These particular directions (along the Euler angles ijj, 9, and <f)) are 
chosen, because the design of the teleoperator is such that it will be 
applying passivating torques to the satellite in these directions. 

These equations can be solved exactly with the help of elliptic 
functions if T = 0. However, in the present analysis, these second- 
order, coupled, differential equations were solved on the digital 
computer using second-order Runge-Kutta method of integration. Various 
despinning and/or detumbling torque functions were tried in each of the 
generalized directions. The torque functions used were: 
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a) zero torque; b) a torque of constant magnitude; c) a torque 
proportional to the angular rate or the coning angle; d) a combination 
of the cases (b) and (c); and e) a triangular torque function. 

The forms of all these torque functions are shown in Figures 5a 
through 5d. 

3.2 Analysis of Motion of the Combined 

Teleoperator-Satellite System 

Under torque-free conditions. In the last section we analyzed the 
motion of the symmetric satellite, moving under the influence of the 
external despinning and detumbling torques applied through the gripping 
arms of the teleoperator. These passivating torques could be applied 
only after the satellite is firmly held in the arms of the teleoperator. 
In most of the previous analyses, it was tacitly assumed that the 
angular momentum vector of the teleoperator coincides exactly with that 
of the satellite. Under these conditions, the passivation of the 
satellite would be a fairly easy task for a symmetric teleoperator- 
satellite combination. However, in practical rendezvous situations, 
it will be most improbable that the teleoperator will approach the 
satellite exactly along its angular momentum vector. In addition, since 
the angular rates of the arms and hand of the teleoperator are to be 
synchronized with the precession and spin rates of the satellite only 
after a short observation pass, there is bound to be some difference in 
the rates and the coning angle of the two bodies. Due to the 
uncertainties involved, it seems appropriate to allow for some deviation 
in the alignment of the angular momentum vector of the teleoperator from 
that of the satellite. In the present analysis, this is achieved by 
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either i) assuming different angular rates for the teleoperator from 
the prescribed precession rate and spin rate of the satellite, or ii) 
different hand angle for the teleoperator from the coning angle of the 
satellite. It is further assumed that the teleoperator mass is an order 
of magnitude greater than the mass of the satellite. Once the tele- 
operator arms grasp the satellite, the satellite starts moving with the 
angular rates prescribed by the arms and its angular momentum vector is 
recalculated. 

With the help of the Lagrangian formulation, the equations of 
motion of the combined teleoperator-satellite mass are derived using 
the expression for the total kinetic energy of the system. This total 
energy is given by 


T_ /. = i[I u^ + I M^ + I 
To/S 2 XX yy z z 
r ^ r 


( 10 ) 


Again, from the geometry of Figure 6, 


u) = ijj- sin 0- sin <|) + 0 cos (^> 

JL M Jl» 

w = sin 0, cos d), - 0, sin d)_ 

y ^1 1 ^1 1 ^1 

0 ) = cos 0- 

z 1 ^1 1 


( 11 ) 

( 12 ) 

(13) 


Substituting these values of o> , o) , and o> into the equation (10) 

X y z 

yields 


1 • *2 

Tto/s 2 ^’^l ®1 ^1 ®1 *^1^ 

1 • *2 

I (^1 sin 0. cos 6_ - 0, sin (b.) 

2 y 1 1 11 1 

' r 

r 


( 14 ) 
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Now, 

body 


utilizing the Euler-Lagrange equations of motion for a rigid 

( 5 ) 


^ ^ _ 9L . 
dt 3 • 

q ^r 


T 

r 


( 15 ) 


as before, and differentiating the expression for the kinetic energy 
of the composite mass with respect to the three generalized 
coordinates, viz., the three Eulerian angles 0^, and will 
yield the equations of motion of the composite teleoperator-satellite 
mass, as follows; 


direction: 


(I sin^ sin^ + I sin^ 0^ cos^ cos^ 0^) 


2*2 * 

+ [21 (sin cos 0^ sin (j>^ 0^ + sin 6^ sin cos 4)^ 


+ 21 (sin 0^ cos 0^ cos^ 0^^ - sin^ 0^ cos sin 4>^) 


-21 (cos 0^ sin 0^ 0^)1 {sin 6^ sin cos <|)^ 0^ 


+ cos 0j^ sin cos <|>^ + 0^^ sin 0^ (cos - sin } 


. 2 


-I {sin sin cos 4>^ 0^ + cos 0^ sin cos 0^' 
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• 2 2 *'i r * *T 

+ 0- sin 0. (cos (|). - sin (Ji.} + I cos 0^ - (}), sin 0, 0,} 

1 A -L liZl 11 11 

r 


= T 




(16) 


direction: 


2 2 ** • 

(I cos $1+1 sin $,) 0. + [I $1 cos 0, sin $- cos $- 
X 1 y 1 1x1 1 1 1 

r j Y r 


- 21 cos $- sin $1 - I cos 0. sin $. cos $, 

x 1 1 1 y ^1 1 ^1 ^1 

r ' r 


+ 21 Sin cos ^^1 ®1 ^1 ^1 


• . • f 2 2 1 • 

+ sin 0j^ {cos - sin $^) - I sin 0^ sin cos 


+ sin 0^ {cos^ - sin^ $^} $j^)] - [I^ sin 0^ cos 0^^ sin^$ 

• • *22 

+ I $, 0, sin $. cos cos 0. + I cos 6_ sin 0, cos 0, 

X 1 1 1 1 1 y ^1 ^1 1 1 

r ' T 


• • * 2 

-I $,01 sin $- cos $, cos0 - I $, sin 0- cos0- 
y 1 1 1 1 z 1 1 1 

'r r 


- '*'1 ®1 """ ®ll = '^ 0 , 

r 1 


( 17 ) 
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~ sin sin cos <{)^ - 0^^ sin cos 4>^ 
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+ 2 i}»j^ 0^ sin 0j { cos^ - sin^ (J>^» 


• 2 2 *2 
+ I ( “ sin 6^ sin cos (j)^ + 6^ sin cos 


- 2 ijj 0 sin 0. {cos^ (J), - sin^ <|).} )] = t 
11 1 1 1 ( 


Since we analyze the motion of the combined teleoperator-satellite 

system in an external torque- free environment, the torques T, , , and 

I'l 9i , 

in the above expressions are identically equal to zero. Also, in 
view of the assumption that the angular rates of the teleoperator are 
imparted instantaneously to the satellite, we do not consider the 
relative motion between the two masses. 


These non-linear, coupled, differential equations were again 
solved on a digital computer using the second order Runge-Kutta method 
of integration. It should be noted at this stage that a preliminary 
design of a single arm teleoperator, used by the Bell Aerospace Co.^^^ 
and shown in the Appendix, is very unstable, even without the spinning 
and processing satellite in its hold. If its arm and the gripping hand 
are spun to the desired angular rates (at which the satellite is 
observed to be processing and spinning), the teleoperator tumbles over 
completely in less than 2 secs as is shown in Chapter IV dealing with 
results. To correct this situation, a modification of the Bell design 
is proposed in this work. Accordingly, the Bell design is modified so 
that the teleoperator consists of four arms and hands, the hands are 
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synchronized to rotate at the same angular rates (although these rates 
may differ from the spin rate of the satellite). The hands rotate in 
such a way that the +z- components of their angular momentum add up and 
the X or y components cancel out each other. This design of the 
teleoperator, as shown in the Appendix, is fotmd to be quite stable 
when its arms and hands are rotating at the prescribed angular rates. 

In all the following results then, only this modified design of the 
teleoperator is considered. 

3.3 Misalignment Between the Angular Momentum Vectors 

As was explained in the previous sections, this misalignment 
between the angular momentum vector of the teleoperator and that of the 
satellite is simulated here by considering the cases .when one or all 
the three variables of the teleoperator differ from those for the 
satellite. The three variables for the teleoperator are: a) the 

angular velocity of the spindle (theoretically matching with the 
precession rate of the satellite) ; b) the angle that the hand makes 
with the +z-axis (matching with the coning angle); and c) the angular 
velocity of the hand (matching with the spin velocity of the satellite) . 

3.4 Initial Conditions 

The misalignment angle is the angle which the +z-axis of the 
composite body (assumed to be coinciding with the +z-axis of the 
symmetric teleoperator) makes with the resultant angular momentum vector 
which remains stationary in space, to calculate this angle and the 
initial location of the body-fixed coordinate frame xyz of the composite 
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body (Fig» 6) with respect to the inertially oriented angular momentum 
vector of the composite body, we assimie that after the satellite is 
gripped firmly in the arms of the teleoperator, it picks up the 
rotational rates imposed on it by the arms. This makes the angular 
momentum vector of the satellite drift from the original +z 
direction by an angle 0^, given by the following relationship. 

I Oi 


9^ = tan 


-1 


xs . xy 
cm ^ 


zs 


0) 


(19) 


Let IJS assume that the hand of the teleoperator may have an initial 
misalignment, then if the actual angle it makes with the +z^axis in 
the satellite-fixed body frame is 0^^ the net misalignment angle which 
the vector makes with the +z-axis in the teleoperator- fixed body 
frame is (Fig. 7a) 


0 


E 
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The resultant momentum vector of the composite body could then 
be easily calculated, and the angle 6^ which the vector makes with 
the + 2 -axis in the teleoperator- fixed body frame, is as shown in 
Figure 7a, given by 


- 0 
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. -1 

sin 


H sin 0„ 
s E 
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r 

Now, if the x-axis of the teleoperator- fixed body frame coincides 
initially with the y-axis of the moving frame, the initial values of the 
Euler angles and angular rates are from Figure 7b as follows. 
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= tt/2 
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(23) 

(24) 

(25) 

(26) 
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Now, we proceed to solve the second order, coupled, differential 
equations of motion from these initial conditions. 
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Chapter IV 
RESULTS 


The results presented in this chapter can be classified in two 
categories. In the first one, only the motion of the satellite under 
the action of the various torque functions is presented. In the 
second case, the motion of the combined teleoperator-satellite mass 
system is presented, both for the case when the angular momentum 
vector of the teleoperator coincides with that of the satellite and 
for the cases when it does not. 

4.1 Motion of a Symmetric Satellite Under 

the Influence of External Passivating Torques 

Figures 5a through 5d show the types of torque functions that 
were applied to the satellite. Although all possible combinations of 
these torque functions, applied in the -direction, 8 -direction, and 
the -direction were tested, only a few relevant results are presented 
in the Figures 8 through 19. To indicate the combination of torque 
functions applied on the three generalized directions, indicators Kl, 
K2, and K3 are noted on each of the Figures 8 through 19. For 
instance, for the case when Kl = 1, K2 = 1, and K3 = 1, it indicates 
that no passivating torques were applied to the satellite; the motion 
is therefore torque-free motion. Similarly, when Kl = 2, K2 = 1, and 
K5 = 2, it indicates that constant torques of magnitudes shown in 
Figure 5a were applied to the satellite in the 4^ and directions, 
whereas no torque was applied in the 0 direction, and so on. 
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The inertia properties of the satellite are described in Table I 
and the initial Euler rates of its motion are tabulated in Table II. 

Table I 

Weight and Inertia Properties of 
Research Applications Module (RAM) Satellite 

Weight = 20,000 lbs. = 9071.84 kg 

Length = 60 ft. = 18.29 m. 

Diameter = 11 ft. = 3.3528 m. 

I = 157,000 lb. ft.^ - 6615.98 kg. 

I = 157,000 lb. ft.^ = 6615.98 kg. m.^ 

yy 

I = 7,050 lb. ft.^ = 297.09 kg. m.^ 

zz 

Table II 

Initial Euler Rates of the RAM Satellite 

ip =0.46 rad/sec. 

(p = 4.0 rad/sec. 

e = 65.9° 

Figure 8 shows the motion of the satellite, when a constant torque 
is applied only in the (f-di recti on. The satellite spin rate is reduced 
to almost zero in 28 secs, and the satellite goes into a flat tumble 
with the coning angle of about 90° and processing at the original rate, 

A satellite is said to be in a state of flat tumble, when it is spinning 
with its angular momentum vector along one of its transverse body axes. 
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Figure 9 shows the motion of the satellite, when a torque 
proportional to the spin rate is applied to the satellite on its spin 
axis. The rate of passivation of the satellite about the spin axis 
is slower, but the coning angle also increases ,at a slower rate than 
in the previous case. Thus, it is seen that the spin velocity of the 
satellite is reduced to a value below 0,5 rad/sec, in about 50 secs., 
whereas the coning angle increases asymptotically to a value of 90^. 

Figure 10 shows the motion when a step-ramp torque is applied only 
about the spin axis. The rate of spin-passivation is faster than in 
the previous case, this being compensated by slightly higher coning 
angles . 

Figure 11 shows the motion when a constant torque is applied only 
on the 0-axis. It is seen that the spin rate is not affected at all; 
the coning angle continues to oscillate, and the precession rate 
though slightly oscillatory in nature keeps increasing over a period 
of time. 

Figure 12 indicates the motion when a constant torque is applied 
only about the ^-axis of the satellite. The spin rate is observed to 
be unaffected, the satellite precession rate is passivated within 
permissible limits within 23 secs., but the satellite assumes a 
position at an angle of about 26° with +z-axis. Figure 13 indicates 
the various angular rates of the satellite when a constant torque acts 
on both the precession and spin axes, with no torque applied on the 
0-axis. The spin rate is passivated within 28 secs, and the precession 
rate within 29 secs. At the end of this period, tho satellite assumes 
a position at about 81° with the +z- axis. Figure 14 shows the motion 
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Figure 9. Motion of Satellite under the Influence of Passivating 
Torques, K1 = 1, K2 « 1, K3 =3. 






Figure 11, Motion of Satellite under the Influence of Passivating 

Torques, K1 “ 1, K2 = 2, K3 « 1, 
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Figure 12. Motion of Satellite under the Influence of Passivating 

Torques, Kl 2, K2 = 1, K3 1. 
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Figure 13. Motion of Satellite under the Influence of Passivating 

Torques, K1 * 2, K2 » 1, K3 = 2* 
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when a constant torque is applied on the ij^-axis and a triangular 
torque function is forced on the spin or ij>-axis. The precession is 
passivated within 24 secs, and the spin rate is also reduced. In 
addition, the coning angle is reduced to about 34^. Figure 15 shows 
the sequence of torque functions applied when it is desired to 
passivate all the three Eulerian variables, 9, and (fi. With the 
sequence shown, the spin is passivated within 28 secs., precession in 
an additional 6 secs., and the deconing is completed in another 4 secs. 
Thus, the satellite is con^letely passivated in 40 secs. It is seen 
from the figure that the decelerations imparted to the satellite are 
fairly regular. Thus, this mode of passivation may be best suited to 
retrieve disabled manned-space capsules. Also, the teleoperator arm 
geometry cannot allow the satellite to reach a pure tumble state. 

Hence, this may also be preferred sequence of passivation for highly 
spinning satellites. Figure 16 shows the motion when a constant torque 
is applied about the precession and spin axes and a step- ramp function 
on the 0-axis. The despinning is completed within 28 secs, and the 
precession and coning is reduced in a total of 40 secs., but by that 
time, a spin velocity of sufficient magnitude in opposite direction is 
imparted to the satellite. So, the passivation is stopped at 28 secs. 
Figure 17 shows the situation when a constant torque is applied on the 
precession axis, a triangular torque fimction on the 0-axis, and no 
torque on the spin axis. The spin rate is totally unaffected and the 
satellite detumbles and decones within 22 secs. The decelerations 
imparted to the satellite, however, are severe. Figure 18 shows the 
motion when a step-ramp torque function is applied on the precession 




0 rad 
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axis, a constant torque on the spin axis, and no torque on the 0-axis. 
The spin is passivated within 28 secs., with the precession rate being 
considerably reduced. The cone angle, however, approaches 90°, thus 
this procedure may not be suitable when the arm geometry does not allow 
the holding of the satellite at 90°. Figure 19 shows the motion when 
a step-ramp torque is applied on the precession axis, a constant 
torque on the 0-axis, and the spin axis is left torque free. The spin 
of the satellite is not affected, and the precession and coning is 
nullified within 23 secs. This procedure may be useful for slowly 
spinning satellites which are processing at high angular rates at 
fairly large coning angles. 

4.2 Motion of the Combined Teleoperator- 

Satellite System in an External Torque- 
Free Environment 

As was pointed out in earlier sections, the asymmetric teleoperator, 

( 2 ") 

as designed by the Bell Aerospace Co.^ was found to tumble over in 
less than 2 secs. Figure 20 shows the angular rates and the coning 
angle of the longitudinal +z-axis. The mass and inertia properties of 
this teleoperator are listed in the Appendix. Since the teleoperator 
should be stable while approaching the satellite, only the improved, 
symmetric design as described in the Appendix is considered. Figure 21 
shows the motion of the teleoperator-satellite combination. In this 
case, the angular velocity of the spindle coincides exactly with the 
precession rate of the satellite, the spin rate of the arm coincides 
with the spin rate of the satellite, and the angle of the arm, 0^, 

■A 

coincides with the coning angle, 0, of the satellite, tt is seen that 
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Figure 18. Motion of Satellite under the Influence of Passivatinj 

Torques, K1 * 4, K2 = 1, K3 = 2. 



Figure 19. Motion of Satellite under the Influence of Passivating 

Torques, K1 = 4, K2 = 2, K3 = 1. 
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Figure 21, Motion of Symmetric Teleoperator After Satellite 

Capture, Case 1. 
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while the precession rate and the spin rate of the combined mass is 
constant, the magnitude of the coning angle 0^ is oscillatory in 
nature. However, the maximum value of 0^ is less than 0.09 radians 
(0.55°), and hence is within the capabilities of the Attitude Control- 
System of the teleoperator Control Unit. Figure 22 shows the variation 
of the coning angle for the combined body, when the spin rate of the 
hand of the teleoperator differs from the spin rate of the satellite 
by 30%, while the precession rate and coning angle are exactly matched. 
It is seen that, again, the precession and spin rates for the combined 
body are constant and the coning angle undergoes irregular oscillations, 
though smaller in magnitude. Figure 23 shows the motion of the 
combined body, when only the angle of the hand differs from the coning 
angle of the satellite by 5°. The oscillations in the magnitude of the 
coning angle 0^ are fairly regular, and the value of the peak 
magnitude decreases slowly with time and is well within the ACS 
capabilities. 

Figure 24 shows the motion of the combined body when the precession 
rate of the spindle differs from the processional velocity of the 
satellite by 10%. Once again, the magnitude of the coning angle 6^ 
for the combined body undergoes oscillations, with peaks slowly 
decreasing with time. The peak value is about 0.53° and well within 
the capacity of the ACS of the teleoperator Control Unit. 

Figure 25 shows the motion of the combined mass when the angular 
rate of the spindle differs by 5%, the angle of hand by 5°, and the 
spin rate of the hand by 30% from the angle and angular rates of the 
satellite. It is seen that the oscillations in the coning angle 0^ of 
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the combined body die out within 25-30 secs., while the angular rates 
remain constant. Figure 26 shows the motion when the spin velocity 
of the spindle differs by 10% from the processional velocity of the 
satellite, the angle of the hand by 5^ from the coning angle of the 
satellite, and the spin velocity of the hand by 30% from the spin 
velocity of the satellite. This is the worst case anticipated during 
the capture. The coning angle 8^ is initially about 0.6^ but again 
dies down within 25-30 secs. 

Thus, it is found that with the modifications indicated in the 
design of the teleoperator, it will be able to capture the satellite 
successfully without itself going into an unstable spin-tumble mode. 
The maximum coning angles are well within the capability of the 
Attitude Control System of the teleoperator Control Unit. 



0 rad 
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Chapter V 
CONCLUSIONS 

A detailed study was presented of the motion of a synnnetric 
satellite under the influence of the passivating torques and of the 
motion of the combined teleoperator -satellite system, with misalignment 
in their angular momentum vectors. It was found that of the many 
torque functions that could be applied, a constant torque would 
passivate the motion in its direction of application faster than other 
forms of torque. For passivating the satellite motion simultaneously 
in all three directions, i.e, if/, 6, and it was found that the best 

method would be to apply a constant torque on all the three axes, 
discontinue the despinning torque when the spin rate goes to zero, 
discontinue the detumbling torque when the precession is nullified, 
and finally decone with the deconing torque. The whole operation would 
take about 40 secs. It is concluded that no general procedure would 
suit all satellites. A particular procedure of passivating would have 
to be developed after abserving the approximate spin and tumble rates of 
the particular satellite to be retrieved. It was also found that with 
modifications in the design of the teleoperator which is mentioned in 
the Appendix, it is possible to successfully capture the disabled 
satellite (before starting the passivation) without the teleoperator going 
into an unstable spin-tumble mode. The maximum coning angle after the 
capture of the satellite, even in the worst case, was of such magnitude 
to be well within the capability of the Attitude Control System of the 
teleoperator Control Unit. 
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APPENDIX 

TELEOPERATOR DESIGN MODIFICATIONS 

The teleoperator and the satellite models used for analyzing 
the motion of the combined satellite -teleoperator system are discussed 
here. The weight and inertia properties of the satellite model used 
and of the teleoperator, as developed in present stages of design, 
are described in Figure 27. However, as is described in Chapter IV, 
this asymmetric teleoperator design appears unstable without complex 
attitude thrusting. As is obvious, there will be a certain interval 
of time from the instant when the teleoperator arm and hand are spun 
up to synchronize with the angular rates of the satellite, to the 
instant when it actually grasps the satellite. Hence, an improvement 
in this design of teleoperator is imperative, and this is illustrated 
in Figure 28. This symmetric teleoperator was found to be completely 
stable when its arm and hand were spun up. The additional inertia 
added to the case or the control unit of this teleoperator was 
necessary in view of the fact that, even though stable by itself, the 
teleoperator (without this additional inertia properties added to the 
case, in the xy plane) would quickly tumble over after holding the 
spinning and processing satellite. This is because of the large asymmetry 
of the combined teleoperator-satellite system in the absence of indicated 
inertia properties. In all the discussions in Chapters III and IV, only 
this modified design of the teleoperator is considered for analyzing the 
motion of the combined teleoperator-satellite system. 




Figure 27. Weight and Inertia Properties of Satellite Despin Model and Teleoperator. 
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Figure 27. (continued) 












